Density-functional studies of spin-orbit splitting in graphene on metals 
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Spin-orbit splitting in graphene on Ni, Au, or Ag (111) substrates was examined on the basis of 
density-functional theory. Graphene grown on the three metals was found to have Rashba splitting 
of a few or several tens of meV. The strong splitting obtained on Au or Ag substrates was mainly 
ascribed to effective hybridization of graphene p z state with Au or Ag d z i states, rather than charge 
transfer as previously proposed. Our work provides theoretical understandings of the metal-induced 
Rashba effect in graphene. 
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Graphene has attracted extensive attention in recent 
years due to its unique and remarkable electronic prop- 
erties, such as gapless-semiconductor band, existence 
of pseudospin, and high electronic mobility at room 
temperature?^— These features are highly desirable for 
the development of next-generation microelectronic and 
spintronic devices 4^ Spin currents in graphene can be 
manipulated using various electronic tactics, in particular 
through exchange and spin-orbit (SO) interactions, which 
are now among the most active research topics in sev- 
eral realms^— The intrinsic SO effect in pure graphene 
layers is nevertheless very weak, 0.1^0.37 meV in flat 
graphene sheet or carbon nanotube* — , insufficient for 
practical use due to the low nuclear charge of the car- 
bon atom. It was hence very exciting when Dedkov et al 
reported an extraordinarily large Rashba 11 SO splitting 
(225 meV) for the ir states of epitaxial graphene layers 
on the Ni (111) substrate through their angle- resolved 
photoemission studies^ It appears that the SO effect or 
hybridization in graphene can be tuned through effective 
electric field across the interface. However, this result 
was challenged by Rader et ati£ who found that the sum 
of Rashba and exchange splitting in the graphene layer on 
either Ni (111) or Co(0001) is less 45 meV. They pointed 
out that the Rashba effect can be strongly enhanced by 
intercalation of one monolayer of Au between graphene 
and Ni (111). Clearly, to discriminate these contradictory 
experimental results and furthermore to understand the 
mechanism of substrate-induced SO splitting in graphene 
are crucial for the progress of graphene physics. 

In this Letter, we report results of density-functional 
theory (DFT)^ calculations for the electronic and mag- 
netic properties of graphene on Ni, Au, or Ag (111) films. 
Interestingly, the Rashba splitting in graphene on Au and 
Ag can be significantly enhanced by strong hybridiza- 
tion of graphene p z state and metal d z i state. The SO 
splitting is found to be almost independent of the charge 
transfer between graphene and its substrates, contradict- 
ing to effective electronic field model proposed in Ref42. 

The electronic structures of graphene on metal (111) 
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FIG. 1: (Color online) (a) Two possible configurations of 
graphene adsorbed on metal (111) substrates. The rhombus 
gives the unit cell along the graphene plane, (b) The energy 
bands of Gr/Ni in the configuration (I) with SO interaction. 
The blue arrow indicates the graphene 7r bands around Y . (c) 
The enlarged 7r bands around F without or with SO interac- 
tions. The exchange and/or SO splitting of the bands are also 
given. The k point in (c) is in the unit of the vector TM. 



substrates, abbreviated as Gr/M (M = Ni, Au, and Ag), 
were calculated by using the VASP code at the level of lo- 
cal spin-density approximation (LSDA)/^ The projector- 
augmented wave (PAW) pseudopotentials were employed 
to describe the effect of core electrons. The equilibrium 
structures were obtained through structural relaxation 
until the Hellmann-Feynman forces were less than 0.05 
eV/A^ The Gr/M systems were modeled by a periodic 
slab geometry, with a vacuum of at least 10 A between 
two neighboring slabs. Each slab contains one graphene 
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layer and N atomic layers of metal. Two adsorption con- 
figurations of graphene on the metal substrates were con- 
sidered, as depicted in Fig. 1(a). We found that configu- 
ration (I) is more stable for graphene on Ni (111), as it 
gives very small lattice mismatch (~1%). For the same 
reason, configuration (II) is preferred for graphene on the 
larger Au and Ag (111) lattice d 17 ' 1 ^ , As a benchmark 
calculation for the SO effect, we first determined the SO 
splitting of the surface states (SS) of the pure Au (111) 
film near the Fermi level (Ep, set as energy zero)J^ Our 
result, AEso ~T00 meV, agrees well with data in RefJ^. 

The band structures of Gr/Ni with N=13^S are given in 
Fig. 1(b) and (c). The direction of Ni magnetization is set 
to be perpendicular to TM. Despite the strong perturba- 
tion from the Ni (111) substrate, one can still easily trace 
several graphene bands in Fig. 1(b), e.g., the graphene a 
and 7r states at -4 and -10 eV in the vicinity of the T point. 
Compared to the band structures of pure graphene, these 
states are spin polarized, and shifted downward in en- 
ergy by about 1.2 and 2.2 eV, respectively. Particularly, 
the feature conical points at K near Ef are destroyed in 
Fig. 1(b), due to broken equivalence of A and B sublat- 
tices through the interaction with Ni. These results are 
in agreement with photoemission measurements^!. 

Now we zoom in to explore spin splitting and SO ef- 
fect of the 7r states of graphene along the -MTM line, 
following the experimental wor k 12 i 13 i 22 . To separate con- 
tributions from different factors, we studied cases cither 
with or without the SO interaction. As illustrated in 
Fig. 1(c), bands without the SO interaction are symmetric 
about the T point and show an induced exchange split- 
ting (Aex) of 30 meV on the magnetic Ni substrate. 
After considering the SO interaction, the energy split- 
ting (Aex+so) contains two parts: exchange and SO 
(Aso). The value of Aso can be extracted through 
Aso — Aex+so — Aex- The linear relationship of 
Aso versus k in Fig. 1(c) indicates that the SO inter- 
action is indeed the Rashba type (Aso = 2auk, where 
an is Rashba strength)] 11 ' 13 The Rashba splitting ob- 
tained from our calculations is about 10 meV, in consis- 
tent with the experimental data in Ref M. In comparison 
with the SO splitting (0.37 meV) in curved graphene 9 , 
the Ni-induced SO splitting in the graphene 7r bands is 
relatively larger. 

The energy bands of Gr/Au in the configuration (II), 
also with a fine lattice match, are given in Fig. 2, where 
different separation (d) between graphene and the Au 
substrate is considered. The 7r bands of graphene are 
accompanied by a gold SS at = -7.3 eV— , labeled as 
'Au SS' in Fig. 2(a). This Au SS is actually localized 
at the vacuum side of the Au slab and hence shows no 
change for different d in Fig. 2(a). The corresponding 
Au SSs at the graphene side move down quickly with 
the decrease of d due to the effect of graphene. At the 
equilibrium geometry (d eg =3.3 A), the SO splitting of 
graphene 7r bands at 0.3 TM is about 21 meV, much 
larger than that on the Ni substrate. For Gr/Au in 
the configuration (II), the A and B sublattices experi- 
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FIG. 2: (Color online) The bands of Gr/Au in the configu- 
ration (II) with N=9 (a) around -7.0 eV along -MTM and 
(b) near Ef along FKM with different d values. The red line 
in (b) is drawn to indicate the Dirac point. The inset is the 
enlarged bands near the Dirac point. 



ence the same environment again. Thus, the graphene 
layer almost restores its unique electronic structure: the 
Dirac cone near Ej?. Charge transfer between graphene 
and the metal and therefore the energy position of the 
Dirac cone can be adjusted through changing d.~ At the 
equilibrium separation, the Dirac point (Fig. 2(b)) is very 
close to E^, indicating almost no charge transfer between 
graphene and Au. When the separation expands/shrinks, 
the Dirac point moves upward/downward, revealing net 
charge transfer between graphene and Au occurs. The 
SO splitting of graphene, however, reduces both ways, 
independent of the enhancement of effective electric field 
in the interface. This contradicts to the effective electric 
field model proposed in Refsi 12 ' 13 for the explanation of 
enhanced SO effect in the systems. The inset in Fig. 2(b) 
shows the SO splitting of the Dirac point. Both the elec- 
tron and hole bands show SO splitting of about 5 meV, 
close to the value (13 meV) given in recent experiment. 24 
The splitting at the Dirac point is induced by the metal, 
as intrinsic SO interaction at this point was predicted to 
be zeroi2£ 

Figure 3 shows the energy bands in large energy range 
for Gr/Au with configuration (II) at equilibrium dis- 
tance. The two sets of parabolas at about -0.2 and -0.4 
eV around the T point are the Au (111) SSs at the two 
sides of the Au slab, with and without graphene, respec- 
tively. They are degenerate at -0.4 eV for the pure Au 
(111) film 19 , but now the SS in the side with the adsorp- 
tion of graphene is pushed up slightly. In comparison 
with Fig. 1(b), the bulk Au 5d bands around T point 
(about -2^-7 eV) are deeper in energy than Ni 3d bands. 
This bestows the probability for Au 5d to strongly inter- 
act with graphene 7r bands, which are in a big interstice 
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FIG. 3: (Color online) The energy bands of Gr/Au(II) with 
SO interaction in large energy range. The blue arrow indicates 
the graphene n bands. The red lines indicate the Dirac point. 



TABLE I: The components of wave functions for graphene 
tv bands at F for the considered systems at equilibrium dis- 
tances. The values are in the scale of the graphene p z state. 
The bold expresses the more stable configuration. The SO 
splitting (in meV) of graphene tt bands in the configuration 
(I) was given at k — 0.25 TM, while in the configuration (II), 
at k — 0.5 TM due to the double lattice vectors. 





Ni (I) 

(N=13) 


Au (II) 

(N=9) 


Au (I) 
(N=12) 


Ag (II) 

(N=9) 


Ag (I) 
(N=12) 


Gr-p z state 


1.0 


1.0 


1.0 


1.0 


1.0 


Ms state 


0.8 


1.3 


0.1 


0.5 


0.7 


M-d z 2 state 


0.0 


2.4 


1.3 


0.0 


3.7 


SO splitting 


6.3 


32.2 


89.0 


1.6 


36.9 



formed mainly by Au 5d states. The situation is very 
different from that in Fig. 1(b), where the graphene tt 
and Ni 3d states are well separated in energy. To un- 
derstand the different SO splittings in Ni and Au cases, 
wave function compositions of the tt bands at T are an- 
alyzed quantitatively. As shown in Table I, the ratio of 
C 2p z :Ni 4s in Gr/Ni (I) « 1.0:0.8; no Ni 3d is involved. 
Since Ni 4s states do not contribute to the SO effect, 
the enhancement of Rashba splitting in Gr/Ni only re- 
sults from the asymmetric potential distribution in the 
two sides of graphene. In contrast, contribution from the 
Au d z 2 state is obvious for Gr/Au (II); the ratio of C 
2p z :Au 6s:Au 5c? is about 1.0:1.3:2.4. Therefore, strong 
hybridization of the metal d z 2 with graphene p z is a key 
factor to produce large SO splitting in graphene tt bands. 
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FIG. 4: (Color online) Left: The bands of Gr/Au (I) along 
-MYM with different d values. Right: The corresponding 
charge density of graphene tt bands at the T point. The den- 
sity is plotted in the plane perpendicular to the interface, 
indicated by the dash line in Fig. 1(a). The upper and lower 
atoms are a carbon and a gold atom, respectively. 



For Gr/Au, configuration (I) is less stable than config- 
uration (II) since the C-C bonds in graphene have to be 
stretched by 17% to match the structure of the Au (111) 
substrate. This stretching gives rise to a much shorter 
equilibrium distance (2.5 A) in the interface, and thus al- 
lows us to understand the effects of adjusting the lateral 
and vertical distances. The whole tt bands of graphene in 
the configuration (I) thus disperse less; and shift upward 
as well, hopefully causing more effective hybridization 
between C 2p z and Au 5d states. Figure 4 gives the SO 
splitting of the graphene tt bands in Gr/Au (I) at d eq 
to be close to 100 meV, 3 times larger than that in Au 
(II) case. In Gr/Au (I), the Au 6s component decreases 
meanwhile the C 2p z and the Au 5d z 2 more effectively hy- 
bridize to each other (see Table I), as expected. Similar 
to the trend in Fig. 2(a) either, increasing or decreasing 
d cause a decrease of the SO splitting in Fig. 4. This 
trend can be rationalized by using the real-space charge 
densities of the graphene tt bands at T in the right pan- 
els of Fig. 4. At the equilibrium distance, a very obvious 
interaction between graphene and metal and asymmetric 
charge distribution above and below the graphene plane 
are observed, corresponding to large SO splitting. When 
d becomes longer or shorter than d eq , the effective mixing 
between C p z and Au d z 2 is weakened. 

Table I also contains the results of Gr/Ag (I) and 
(II). Again, the SO splitting in the configuration (II) is 
much less than that in (I) due to weakened hybridization. 
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TABLE II: The SO splitting of n bands of graphene on differ- 
ent metal substrates with N=l and 6, respectively. The value 
is given at k — 0.25 or 0.5 TM, as stated in Table I. The bold 
indicates the more stable configuration. 



meV 


Ni (I) 


Au (II) 


Au (I) 


Ag (II) 


Ag (I) 


N = 1 


6.7 


11.2 


35.4 


2.6 


6.7 


N = 6 


6.3 


33.4 


88.0 


2.4 


41.8 



While the graphene p z state has strong interaction with 
the Ag d z 2 state in the configuration (I) , only Ag s state 
is involved in the more stable Gr/Ag (II), similar to Ni. 
Therefore, heavy metals may not always produce large 
SO splitting in graphene. 

Finally, we also explored the effect of the thickness of 
metal films. The SO splitting of graphene ir bands for 
Gr/Au (I) and (II) and Gr/Ag (I) with N=6 in Table 
II are almost the same as the values listed in Table I, 
respectively. For the rest two cases: Gr/Ni (I) and Gr/Ag 
(II), only one monolayer of metal substrate is enough to 
give a saturated SO splitting of graphene. Since Au SSs 
usually extend several atomic layers into the bulk 23 , a few 
layers of metals are needed to obtain the saturated SO 
splitting. Nevertheless, the SO splitting for graphene on 



Au mono- and bi-layer films is already large, explaining 
why one Au atomic layer intercalated between graphene 
and Ni (111) can cause a substantial Rashba effect in the 
experiment^ 

In conclusion, we investigated what determines the SO 
splitting of graphene bands on Ni, Au, or Ag (111) sub- 
strates through first-principles calculations. While the 
Rashba splitting for Gr/Au is sizeable, the effect of Ni 
is very limited. The hybridization between graphene p z 
and metal d z i states is identified as the chief factor for 
the enhancement of the SO effect in graphene. This re- 
quires not only large SO strength from metal atoms, but 
also effective overlap of metal d and graphene states in 
energy. A few atomic layers of metals are sufficient to 
produce saturated strong Rashba splitting in graphene. 
Our findings point out a direction for the manipulation 
of SO strength in graphene that is needed for the devel- 
opment of spintronic materials and devices. 
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